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Abstract 

We calculate the full asymptotic expansion of boundary blow-up so- 
lutions (see equation JT} below), for any nonlinearity /. Our approach 
enables us to state sharp qualitative results regarding uniqueness and ra- 
dial symmetry of solutions, as well as a characterization of nonlinearities 
for which the blow-up rate is universal. Lastly, we study in more detail 
the standard nonlinearities /(«) = u p , p > 1. 

1 Introduction 

Let B denote the unit ball of R N , N > 1 and let / G C(R). We study the 
equation 



(1) 



Au = f(u) in B, 

u — +oo on dB, 



where the boundary condition is understood in the sense that 
lim u(x) = +oo for all xo € dB 

and where / is assumed to be positive at infinity, in the sense that 

(2) 3 a <E R s.t. f(a) > and f(t) > for t > a. 

A function u satisfying JTJ is called a boundary blow-up solution or simply 
a large solution. Existence of a solution of (JTJ is equivalent to the so-called 
Keller-Osserman condition : 



r+°° dt r* 

(3) J 7f ^< + oo, where F{ t)= ^ 



f(s) ds. 



1 
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For a proof of this fact, see the seminal works of J.B. Keller [6] and R. 
Osserman [8] for the case of monotone /, as well as [4] for the general 
case. From here on, we always assume that Q holds. 

Our goal here is to study asymptotics, uniqueness and symmetry prop- 
erties of solutions. Our approach improves known results in at least two 
directions : firstly, aside the necessary condition ((3]), we need not make any 
additional assumption on / to obtain the sharp asymptotics of solutions. 
Secondly, we obtain the complete asymptotic expansion of solutions, to 
all orders. Here is a summary of our findings. 

Theorem 1.1 Let f G C(R) and assume ((2|, ([3]) hold. Consider two 
solutions ui , 1*2 of {TJ . Then, 

lim ui(x) — U2(x) = 0, for all xo G dB. 

x — >XQ ,x£B 

More precisely, there exists a constant C — C(ui, U2, N, F) > 0, such that 
for all x G B, 

f +oc dt 

(4) \ Ul (x)-u,(x)\<C —dt. 

Ju 2 (x) * W 

In addition, 

(5) |.F(ui) - F(u 2 )\ G L°°(B). 

Estimates on the gradient of solutions can be obtained for a restricted 
class of nonlinearities, namely 

Theorem 1.2 Let f G C(R) and assume © hold. Assume in ad- 
dition that f is increasing up to a linear perturbation i. e. there exists an 
increasing function f and a constant K such that 

(6) /(*) = f(t) - Kt, 

for all t£l. Consider two solutions Mi,M2 of HJ. Then, 

(7) |V(u 1 -u 2 )|GL°°(B). 

As a direct consequence of Theorem 11.11 we obtain the following unique- 
ness result. 

Corollary 1.3 Let f G C(R) and assume © hold. Assume in addi- 
tion that f is nondecreasing. Then, there exists a unique large solution of 
©■ 

Remark 1.4 Many uniqueness theorems have been established in the lit- 
erature (see e.g. the survey [T]^, and they hold for a general class of 
bounded domains O,. However, in all of these results, additional assump- 
tions on f are needed, such as convexity. 

Proof of Corollary 11.31 . Let Ui,U2 denote two large solutions. It 
suffices to prove that u\ < U2- Assume this is not the case and let uj = 
{x G B : w(x) > 0} 0, where w = ui — u-z. Working if necessary on 
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a connected component of w, we may always assume that to is connected. 
Using Theorem ll.il we see that w solves the equation 

j Aw = f(ui) - f(u 2 ) > in u, 

w = on 9u>. 

By the Maximum Principle, w < in u, a contradiction. □ 

When / is not increasing, uniqueness fails in general. One may ask 
however whether all solutions of fl]) are radial. H. Brezis made the fol- 
lowing conjecture. 

Conjecture 1.5 ([3]) Let f G C 1 (R) denote a function such that 
([3j) /ioZd. Then, every solution of (JXJ) is radially symmetric. 

A. Porretta and L. Veron proved that this is indeed the case under the ad- 
ditional assumption that / is asymptotically convex (see [9])- We improve 
their result as follows. 

Corollary 1.6 Let f G C 1 (R) and assume l|2"]l. © hold. Let u denote a 
solution of JTJ). Assume in addition that, up to a linear perturbation, f 
is increasing (i.e. ([6]) holds for some nondecreasing function f and some 
constant K). Then, u is radially symmetric. Furthermore, ^ > in 
B\{0}. 

Remark 1.7 In the setting of the classical symmetry result of B. Gidas, 
W. M. Ni and L. Nirenberg (see ((6| is also assumed in order to 

prove symmetry. In the same article, the authors give an example of 
a nonlinearity f failing © for which there do exist nonradial solutions 
of the equation. In the context of large solutions, we do not have such 
a counter-example. In fact, we expect that none exists i.e. we believe 
that Conjecture \1.5\ holds true. But at this stage, we do not even know 
whether radial symmetry remains true for simple nonlinearities such as 
f(u) = u 2 (l + sinw). 

Corollary 11.61 is a direct consequence of the moving plane method and 
Theorem O: 

Proof of Corollary 11.61 . Let U denote a radial solution of |T}. It 
follows from (fSJl that U is a nondecreasing function of r = \x\ for r close 
to 1~ and ^r(r) — > +oo as r — > 1~ . By 1(7}. we conclude that any solution 
u of fTJ satisfies §7(2;) — > +00 as x — > dB, while the tangential part of 
the gradient of u remains bounded. We then apply Theorem 2.1 in [9]. □ 

In addition to the relative asymptotic information given by )4|, |5) 
and 0, the exact asymptotic expansion of a solution can be calculated 
to all orders, as follows. 

Theorem 1.8 Let f G C(R) and assume ©, © hold. Let Uo G K, 
/ = [Uo, +00) and let vo be the function defined for u G / by 

(8) v {u) = y/2F(u). 

Consider the Banach space 

X = {v G C(7;K) : 3M > such that \v\ < Mv }, 



cdl9mars2009.tex 



March 20, 2009 



4 



endowed with the norm \\v\\ — supj |u/i>o|. If the constant Uo is chosen 
sufficiently large, there exists p £ (0, 1) such that there exists a unique 
solution v £ B(vq, p) C X of the integral equation 



(9) 



v{u) = < 2 IF (u) - (N - 1) 



- dt 

r 



we I, 



where r = r(u, v) is given for u 6 I, t) £ B(vo, p) by 



(10) 



r(u, v) = 1 



Ju V 



(It. 



In addition, v is the limit in X of (vk) defined for k = by (JSJ) and for 
k > 1, by 



(11) v h (u) = 



2 F(u) - (AT - 1) 



Ufc-l 



r.'o 



^ _ r+oo 1 



and t/ie sequence vu is asymptotic to v i.e. as u — > +oo, 
Wfc+i(«) = «k(u) + o(v k (u)) 

and given any k £ N, 

l>(w) = t)fe('ii) + 0(vk+i(u) - Vk{u)). 

Let now u denote any solution of JTJ) and fix ro £ (0, 1) such that u(x) > 
Uo for \x\ > ro- For k > 0, define Uk for r > ro as t/ie unique solution of 



(12) 



dw fe 

-3— = V k (life ) 

ar 

lim itfc(r) = +00, 

^ r— 



where Vk is given by (|11|) , Then, as r — > 1 , 

r k + i{r) dM _ ( r +o ° du 

Ju k (r) VO \Ju k (r) «D 

and given any k £ N, u;e /lane as a: — -> 9-B, 



(13) 



u(a;) du 



Theorem 1 1 .8 1 enables one to calculate (implicitely) the asymptotic expan- 
sion of a solution term by term. But how many terms in this expansion 
are singular? This is what we discuss in our last set of results. 

We begin with the simplest class of nonlinearities /, those for which 
only one term in the expansion is singular, namely the function uo defined 
by and (|12|). It turns out, as A.C. Lazer and P.J. Mc Kenna first 
demonstrated (see [7]), that in this case «o(l — d(x)) is the only singular 
term in the asymptotics of any blow-up solution on any smoothly bounded 
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domain Q, C and for any dimension N > 1, where d(x) denotes the 
distance of a point x £ Q. to the boundary of Q,. In other words, the 
blow-up rate is universal. The question is now to determine for which 
nonlinearities /, this universal blow-up occurs. We characterize these 
nonlinearities as follows: 

Theorem 1.9 Let C denote a bounded domain satisfying an inner 
and an outer sphere condition at each point of its boundary. Let f £ C(M), 
assume ([2J, Q hold and consider the equation 



(14) 

Assume 
(15) 



Ait = f(u) in Q, 

u = +oo on dQ.. 



r+°° fcy/2Fds 



Then, any solution of (|14[1 satisfies 

(16) lim u(x) - m (1 - d(x)) = 0, 

x— >dQ 

where d(x) — dist(x, <9f2) and uo is defined by |[BJ, l|12|). 

We also have the following partial converse statement : if 

r+oo r t J2F ds 

(it) i^ 1 Vmu)j u ^^ r «ft>o, 

then (|16|) always fails. 

Remark 1.10 To our knowledge, (|15[) improves upon all known condi- 
tions for i|16[) to hold (see in particular [7] and [2],) . Despite its unap- 
pealing technical looks, (|15p uses information on the asymptotics of F only 
(in particular, no direct information on f is required). Nonlinearities such 
that F(u) ~ e u or F(u) ~ u p , p > 4 as m — > +oo qualify. For F(u) ~ u 4 , 
(|17l) /to£ds and so t/ie conclusion (| lCjp /ai£s. 

Remark 1.11 Condition (1151) can 6e weakened to : 

rt 



lim x/2F(m ) / 

f—* 1 inn 



J ° „Y — 



where uo = uo(r) is defined by (|12l) . Similarly, (|17[) con 6e weakened to 

/+°o r* A /2F ds 
( 2F)3/2 * > - 

As an immediate corollary, we obtain uniqueness on general domains, 
whenever only one singular term appears: 

Corollary 1.12 Assume (|15|) . // in addition, f is nondecr easing, then 
the solution of (|14|) is unique. 
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Proof. Simply repeat the proof of Corollary 11.3 



□ 



More than one term can be present in the asymptotic expansion of u. 
Finding all the (singular) terms in this expansion is of staggering algebraic 
complexity. To illustrate this, we provide the first three terms (in implicit 
form) . 

Proposition 1.13 Let ui be defined by (|12[1 for k — 2. Let also Ri, R2, 
7?3 denote three real-valued functions defined for U £ R sufficiently large 
by 



Ro(U) = 
( 



Ju 



du 

72F' 



(AT- 



Ri(U) = (N — 1) 

R 2 (U) = (N — 1) x 
fV2F ds 



f 

Ju 



J u V2Fdt 

(2F)3/ 2 



du, 



<2F 



+ V2F 



ds 

J2F 



dt- 



V 



hn-1) {rv^Fdtj 

4 2F 



(2F)3/ 2 



Then, for all r £ (0, 1), r close to 1, there holds 

1 - r = R (u 2 (r)) + Ri(u 2 (r)) + R 2 (u 2 (r))(l + o(l)). 

For specific nonlinearities, it is possible to invert the above identity. 
This is what we do for f(u) — u p , p > 1: 

Proposition 1.14 Let p > 1 and f(u) = u p , foru > 0. Then, the unique 
positive solution of {TJ satisfies as r — > 1~ , 

[2/<p-l)l 

u = cTp^ t 22 a k d k + o(l), 

fc=0 

where d(r) = 1 — r for r £ (0, 1), and where each at G R depends on N 
and p only. 

Remark 1.15 Note that the above result remains true for any nonlinear- 



ity f such that, for some positive constant c, F(u) ~ cu 
(and any solution of the equation). 



p+i 



+00 



Outline of the paper 

1. In the next section, we show that any solution u of Q can be 
squeezed between two radial solutions U and V i.e. the inequal- 
ity U < u < V holds throughout B. 

2. Thanks to this result, we need only find the asymptotics of radial 
solutions to prove Theorem ll.il This is what we do in Section [3] 

3. To obtain gradient estimates, the squeezing technique is insufficient 
and more work is needed. In Sectionf?] we estimate tangential deriva- 
tives via a standard comparison argument, while we gain control over 
the radial component through a more delicate potential theoretic ar- 
gument. 
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4. Section [5] is dedicated to the proof of Theorem 11.81 that is we estab- 
lish an algorithm for computing the asymptotics of solutions to all 
orders. 

5. In Section [6] we characterize nonlinearities for which the blow-up 
rate is universal. 

6. At last, Sections [7] and [5] contain the tedious calculations of the first 
three terms of the asymptotic expansion of u in implicit form for 
general /, and of all terms explicitely for f(u) — v? . 

Notation 

Throughout this paper, the letter C denotes a generic constant, the value 
of which is immaterial. In the last section of the paper, we use the symbol 
Cfc to denote a quantity indexed by an integer k, thought of being "constant 
for fixed k" , the value of which is again immaterial. 



2 Ordering solutions 

In this section, we prove that any solution of the equation is bounded 
above and below by radial blow-up solutions. To do so, we impose the 
following additional condition: g(t) := f(—t) satisfies ([2]) and 

(18) J -^g==+oo, where G'(t) = f(-t). 



Remark 2.1 Note that ()18[) is not restrictive. Indeed, if u denotes a 
solution of |TJ and m = rnins u, then u also solves UJ with nonlinearity 
f defined for u € R by 



f(m) + (m — u) if u <m 
f(u) ifu>m. 



Then, f clearly satisfies (|18|) . 

We now proceed through a series of three lemmas. 

Lemma 2.2 Assume (|18[) holds. For M G R sufficiently large, there 
exists a radial function v £ C 2 {B) D C(B) satisfying 

Av > f(v) in B 

and such that 

v < -M in B. 

Proof. Let g(t) = f(—t) for t G R and let a > M a parameter to be fixed 
later on. Since g satisfies ((2j), we may always assume that g(t) > for 
t > M. Let now w denote a solution of 



(19) 
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Claim. There exists a > M sufficently large, such that w(l) > M. 

Note that w is nonincreasing in the set {t : w(t) > M}. We distin- 
guish two cases. 
Case 1. w > M. 

In this case, w is defined on all of R + . In particular, w(l) > M, as 
desired. 

Case 2. There exists R > such that u>(i?) = M. 

In this case, since w is nonincreasing in (0,R), we just need to prove 
that R > 1. To do so, multiply (|19(l by — iw and integrate between and 
r G (0, R): 

-w' = y/2(G(a)-G(w)), 
where G is an antiderivative of g. Integrate again between and R: 

dt f - W dr = R. 



m v/2(G(a) - G{t)) Jo s/2{G{a) - G(w)) 
Now, 

/•« f lf rG- 1 (G(a)/2) 

R= => / => 

Jm v/2(G(a) - G(t)) 7m \/2(G(a) - G(f)) 

G- 1 (G(a)/2) rft 



> 



v/2G(t) 



By (|18|) . we deduce that _R > 1 for a sufficently large. We have just proved 
that w\, !j > w(l) > M and the claim is proved. 

It follows that the function « defined for x G B by = — w(\x\), is 
the desired subsolution. □ 

Lemma 2.3 Let / G G(R) and assume Q Zioirf. Asjumeu G G(B) 
satisfies 

Av>f(v) mB. 
Then, there exists a radial large solution V of {TJ such that V > v. 

Proof. Let v := N. Then, v and v are respectively a sub and supersolu- 
tion of 



(20) 



Av = /(«) in B, 

v — N on 9£?, 



provided N is chosen so large that N > \\v\\l°^(b) and f(N) > 0. Futher- 
more, < u in S for such values of AT. By the method of sub and 
supersolutions (see e.g. Proposition 2.1 in [4]), there exists a minimal 
solution Vn of (|20[) such that iV > Vn > «■ Note that Vjv is radial, as 
follows from the classical symmetry result of Gidas, Ni and Nirenberg (see 
[5]) when / is locally Lipschitz, or simply from the fact that Vn is mini- 
mal, hence radial, for / merely continuous. Also, since Vjv is minimal, we 
have that the sequence (Vjv) is nondecreasing with respect to (apply 
e.g. the Minimality Principle, Corollary 2.2, in [4]). 

It turns out that the sequence (Vjv) is uniformly bounded on compact 
sets of B. Indeed, fix Ri < 1. There exists a solution U blowing up on the 
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boundary of the ball of radius 1 and satisfying U > v_ in Br 1 , see Remark 
2.9 in [J. By minimality, v < Vn < U in Br 1 , whence (Vn) is uniformly 
bounded on Br 2 , for any given R2 < Ri . 

We have just proved that each Vn is radial and that the sequence (Vn) 
is nondecreasing and bounded on compact subsets of B. By standard 
elliptic regularity, it follows that (Vn) converges to a radial solution V of 
HJ, such that V > v in B. □ 

Lemma 2.4 Assume ((3} and (|18[) hold. Let u denote a solution of ifl]). 
Then, there exist two radial functions U, V solving {TJ such that 

U <u<V in B. 

Proof. Let —M denote the minimum value of u and let v denote the 
subsolution given by Lemma 12.21 In particular, v < u. By Lemma 12.31 
there exists a solution U > v of ([IJ and we may asssume that U is the 
minimal solution relative to v i.e. given any other solution u > v of (fTJ), 
U < u. In particular, [/ < u. It remains to construct a radial solution V" 
of |T|) such that u < V. To do so, we fix R < 1. By Lemma 12.31 letting 
v = u\ B , there exists a radial solution v = Vr of 




= f(v) in 
d = -j-00 on 8Br, 



such that Vr > u in Since Vr is constructed as the monotone limit of 
minimal solutions Vn (see the proof of the previous lemma), one can easily 
check that the mapping R 1— > Vr is nonincreasing (hence automatically 
bounded on compact sets of B). Hence, as R — > 1, Vr converges to a 
solution V of ([T]), which is radial and satisfies V > u in B, as desired. □ 



3 Asymptotics of radial solutions 

Our next result establishes that the asymptotic expansion of a radial blow- 
up solution is unique. More precisely, consider the one-dimensional prob- 
lem 

(22) = /(£), r <l, *(r)^+»asr^r. 

All solutions are given implicitely by 

r+°° d s 

I — , = 1 — r, where F = f . 

h V2F(s) 

We recall the following fact, first observed by C. Bandle and M. Marcus 
in [2] : 

Remark 3.1 Let tf> and cf) c denote two solutions of (|22[) corresponding to 
the antiderivatives F and F + c, respectively. Then 4>(r) — 4>c(r) — > as 
r -> 1" . 

We improve this result in the following way. 
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Theorem 3.2 Let N > 1 and let «i,Ua denote two strictly increasing 
functions solving 



(23) 
T/ien, 



d 2 u N -ldu „, , 

TT + T" = /(")' r<1 > 

dr 2 r dr 

lim w(r) = +oo. 



|wi(r)-w 2 (r)| < C* / dt 



dt 

u. 2 (r) F ( t ) 



In addition, the quantity \F(ui) — _F(m2)| is bounded. 

Remark 3.3 Clearly, Theorem ] 1 . 1\ follows as a direct consequence of Re- 
mark [2A\ Lemma\2.4\ and Theorem \3.1& 



Proof. We want to think of the second term on the left-hand side of 
equation ([23j) as a lower order perturbation as r — * 1. So, we integrate 
(|23|) in the same way we would solve Q22| l. namely we let v = du/dr and 
multiply the equation by v. We get 

d fv 2 \ N - 1 2 d „ 

— — H w = — (F(u)) . 

dr V 2 y r dr v w; 

We define the resulting error term by 

(24) g:=~+F(u), 

which, seen as function of r, satisfies the differential equation 

(25) d -f = ^lv\ 

dr r 

Since u is a strictly increasing function, the change of independent variable 

u — u(r) is valid. Thinking of g as a function of the variable it, we have 
dg = dst = ids and gQ 

au ar au v ar 

(26) £ = ^l v . 

au r 



Since (|23[) holds for r close to 1, the above equation holds for u in a 
neighborhood of +oo. Solving (124[) for v, we finally obtain 



TV — 1 N 



(27) 



v= ^2{F{u)-g), 

i r 

u =l/v = (2(F(u)-g))- 1 / 2 . 



We start by calculating the leading asymptotic behaviour of g at +oo : 
Lemma 3.4 

lim 4t=0. 

In addition, 

(28) lim — ^^ U )„, . — 1, where G is any antiderivative of V2F. 
y ' u-f+oo (N - l)G(u) 
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Proof. First, we claim that 

(29) lim #4=0- 

u^+oo F[U) 



Indeed, fix e > and recalling that ((3]) holds, choose M > so large that 

r + oc 
J At 

that 



f„ + r °° ; dt < e. By the definition of G, there exists a constant Cm such 

JM % /2F(t) 



G(u) =C M + / y/2F(t)dt. 

JM 

Since F is nondecreasing it follows that 

G{u) < Cm + 2F(u) / < Cm + 2eF(w) 

./m \/2F(t) 



Dividing by F(u) and letting u — > +oo, (|29l) follows. Next, we claim that 
(30) lim |^-=0. 

Note that by (|26[) . <?(m) is increasing, thus bounded below by a constant 
c as m — + +oo. Hence, by ([27)1 , 

N — 1 



^ < —— -V2(F&) ~c)< 2{N - l)vW - c), 
cm r 

where the last inequality holds if r > 1/2 i.e. if it is sufficiently large. 
Integrating on a given interval (uq,u), we obtain 



< 9{u) < g(uo) + 2(N - 1) / y/2(F(t) - c) dt 



Using (1291) and the fact that limt_,+ 00 \/ 2F ( t '>_ _ i we deduce (1301) . 

' — ' V 2 ( F (*)- c ) 

Now that 1)301) has been established, we return to (|27[) and infer that 
given e > 0, we have for sufficiently large u, 

^ > ^-V2(l-e)F(u)> (N - 1W2(1 - e)F(u) 
du r 

and 

g < ^-lv^+e)F(u) < ^-l V 2(l + e)F(u). 
Integrating the above, we finally obtain for large u, 

(l-e)(JV-l) / ^2Fif)dt <g{u)-g{u ) < (l+e) 3/2 (N-l) / ^2F\F) 

JUQ J UQ 

and (|28[) follows. The fact that g(u) — > +oo as u — > +oo follows automat- 
ically. □ 

Next, we prove that given two solutions ui, U2, the corresponding error 
terms gi , g2 given by (|24[) differ by a bounded quantity. 



dt 



cdl9mars2009.tex 



March 20, 2009 



12 



Lemma 3.5 Let u\,u 2 denote two solutions of (|23l) . Introduce m = 
and ^ 

3i = -y + -FXO, /or i = 1,2. 

Then, gi — 172 is bounded. 

Proof. We have seen in the course of the proof of Lemma 13.41 that each 
gi (seen as a function of a variable u lying in some interval (Mi, +00)) 
solves the differential equation (|27[) , Letting 10 = pi — 32, it follows that 

dw N - 1 



(^2(fM -ffi)-N/2(f>)-<72)) 



du r 

2(JV-1) 1 



r y/2(F(u) - 51) + v/2(F(«) - g 2 ) 
So, for some constants C, uo, 

/ f u 2(N-l) 1 
w = C exp - / -i i ■ d£ 

V Jua r y/2(F(t) - gi ) + y/2(F{t) - g 2 ) 
whence \w\ < \C\ for u > uo. □ 

Completion of the proof of Theorem 13.21 

Let ui,u 2 denote two solutions of (|23p . By (|27|l . each Uj, i = 1,2, 
solves 

dtti / dr 

Integrating, we obtain 

dt = 1 - r = / dt. 



= 1. 



J U 



«, ^2(F(t) - ah) i U2 v/2(F(i) - 02 

Without loss of generality, for a given r we may assume 112(f) > tii(r). 
Split the left-hand side integral : f + °° = f" 2 + f + °°. It follows that 



p i -2 



dt= I - dt 



y/2(F(t)- gi ) Ju 2 \y/2(F(t) - g 2 ) y/2(F(t)- 9l ) 

+°° y/2(F(t) - gi ) - y/2(F(t) - g 2 ) it 
2 y/2(F{t)- gi )y/2(F(t)-g 2 ) 

g ' 2 - 91 dt 



'u 2 y/2(F(t) - gi)^2(F(t) - g 2 ) (y/2(F(t) - gi) + ^2{F{t) - g 2 )) 

Recall that by Lemma [3.4l gi = o(F) as t — > +00. Recall also that g 2 — gi 
is bounded. So, for sufficiently large values of u 2 , we deduce 



f u 2 1 /H 

(31) J -7wm dt - c 



dt 



F(t) 3 /* 
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Since F is increasing, it follows that 

dt 



o< U2 ~ Ul 



/•"2 1 r-t 

<C -== dt<C 
J «, \/2F(tj Juo 



V^M~ J U1 y/Wtfj - L 2 F(t) 3/2 

C f + °° dt 



< 



v/fR Ju 2 F(t) ■ 

Hence, 

f +OQ dt 

as stated in Theorem 13.21 It remains to prove <j5j - Without loss of gener- 
ality, we assume ui(r) < u%(r) so 

r 2 dt _ r + °° dt r +o ° dt 

r +cc dt r +oc dt 

■J u 



y/F{t - (U2 ~ til)) Ju 2 y/Fj£) 

L 2 ^F(t)F{t - {u 2 - iti)) 



yfF{t)F{t - (u 2 - Ml)) + - ("2 - Ml))) 



f + °° dt 
> (F(u 2 ) - F( U1 )) / 



F(t)3/2" 

Recalling ([3l]l. © follows. □ 



4 Gradient estimates 

Proof of Theorem 11.21 Let w = ui — U2 denote the difference of 
two solutions. Without loss of generality, we may assume that U2 is the 
minimal solution of {l}, so that ui > 112 and U2 is radial. 
Step 1 : estimate of tangential derivatives 

We begin by proving that any tangential derivative of w is bounded. 
Since the problem is invariant under rotation and since 112 is radial, we 
need only show that (r, 0, . . . , 0) remains bounded as r — > 1~ . Given 
x = (xi, X2, x') £ B and 6 > small, we denote by xe = (^icosS — 
X2 sin 8, x\ sin # + X2 cos a;') the image of x under the rotation of angle 6 
above the xi-axis in the (x\,X2) plane. By the rotation invariance of the 
Laplace operator, the function ug defined for x £ B by u$(x) = ui(xe), 
solves (flj. Using @ and assumption ©, we deduce that Wg = ui — ug 
solves 



Awg + Kwg = /(tii) — f(ug) in B, 

wg — on dB. 
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By the Maximum Principle on small domains, there exists Ro £ (0, 1) such 
that the operator L — A + K is coercive on B \ Br q . As a consequence, 
we claim that there exists a constant C > such that for all x £ B \ Br , 



(33) 



\wg(x)\ < C sup \wg\ 



OB 



Ho 



Let indeed ( > denote the solution of 
AC + KC, = 

C = i 
C = o 



in B \ B Ro 
on dB Ro 
on dB. 



±sup eSfi |?f£i|C are respectively non- 



We shall prove that z := 
positive and nonnegative, which implies that (|33p holds for the constant 
C = HCIloo- We work say with z + and assume by contradiction that the 
open set u = {x £ B \ Br : z + (x) > 0} is non-empty. Restricting the 
analysis to a connected component, we have 



Az+ + Kz + = /(mi) - f(u g ) > incj 
z + < on duj. 



By the Maximum Principle, we conclude that z + < in to, a contradiction. 
We have thus proved (|33|1 , Since iti £ C 1 (Bh ), we deduce that for some 
constant C > and all a; £ B \ S flo , 



\w e (x)\ < C6. 



Applying the above inequality at the point x 
and letting 6 — > 0, we finally deduce that 



•0) 



< c 



(r,0,...,0), r £ (i? ,l) 
for all r£(R ,l), 



as desired. 

Step 2 : estimate of the radial derivative 

It remains to control dw/dr. Fix R £ (0, 1). Let Gn(x,y) denote Green's 
function in the ball of radius R. Then, for x £ Br, 



(34) 



w(x) 



dG 



R 



I dB R 9v y 

=: W!(x) 



(x,-)wda + GR(x,-)(f(ui) - f(u 2 )) dy 

+ W2(x). 



We want to let R — * 1 in the above identity. To do so, we first observe 
that wi is harmonic. By the Maximum Principle, |wi| < HujU^oojgB^). 
By estimate Q, we conclude that wi — » as R — * 1. To estimate u>2, we 
need the following crucial estimate : 

Lemma 4.1 Assume (El). Then, 



\f{ui)-f(u 2 )\(r,0)dr< 



sup 

9 g S N-l J 
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We shall also need the following elementary estimates. 

Lemma 4.2 There exists a constant C > such that for all 1/2 < r, R < 
1 and all x, y £ Br, 



(35) 



G R ( X ,y)=R"- N G 1 (|,|) 



(36) 



Gi(av) da < 1 



dB r 



^(X •) 

a\x\ 



da <C 



We postpone the proofs of the above two lemmas and return to Q34 
Using polar coordinates, 

w 2 (x) = G R (x,-)(f(ui) - f(u 2 )) dy 



dB T 



G R (x,-)(f(ui)- f(u 2 ))da) dr 



By Lemmas 14.11 and 14.21 we may easily pass to the limit in the above 
expression as R — > 1, so 

w{x)= f Gi(x,-)(/(ui)- f(u 2 ))dy 



Using again Lemmas 14. II and 14.21 we also have that w is differentiate in 
the r = | a; | variable and 



Using polar coordinates again and Lemmas 14. II and 14.21 we finally obtain 
dw 



< C+ 



dr 



+ sup 

r€(l/2,l) \JdB r 



8Gi 



(*,•) 



da ) sup 



i (/ ICf( ul )-/MIM) dr 



< c. 



ft only remains to prove Lemmas 14.11 and 14.21 

Proof of Lemma 14. II We first deal with the case where ui, u 2 are radial 
and Mi > u 2 . By assumption ((6]), we have 

|/(ui)-/(ua)| dr <f o - /(«a)) dr + K\\ Ul -u 2 \\ L oo (B) . 

Using Q, we see that u\ — u 2 is bounded and so it remains to estimate 
/ (ui) — f(u 2 ). By (f27|) . each iti, i = 1, 2, solves 



^(^(uO-sO 
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We also know by Lemma 13.41 that gt = o(F(ui)). So, 

.. dui/dr 

lim — ■ = 1. 

Using this fact, as well as Lemma 13.51 and (|3J, we obtain for R £ (1/2, 1), 

f R (f(ui)-f(u 2 ))dr< f R (/(«i) - /(«a)) dr + iri|m-M2|U~ (s) 

< C A/K) - /(^)) dr + C 

Jo v2F(ui) 

\" ( R !': )( dU2,dr dui / dr ) dr + c 

< C (yF(vi) - VF(u 2 )) (R) + C+ 

C r /M ( W£±J* _ y/2(FK)- gl )\ ^ 
Jo '\ sfW(u2} y/Wfc) ) 

<c{ lM-^= ) (R ) + C+ 

, n f R fl ^ y/F(M 2 )F(m) - g 2 F(m) - ^/F( Ul )F(u 2 ) - 9l F(u 2 ) 

+ 6 / }{u 2 )— = ar 

Jo v^W%0 

||F(m)-F(« 2 )|U~ (S )+C+ 



+ 



, f R (M giF(u 2 ) - g 2 F( Ul ) ^ 



o y/F(u 3 )F(ui) y/F(u 2 )F( Ul ) 
<C + C f R f(U2) (9i-g2)F(u 2 )+g 2 (F(u 2 )-F( Ul )) ^ 
Jo ^F{u 2 )F{ Ul ) v /F(u 2 )F(«i) 

<c + cn ffl - 52 || L oo (fl) / fl iM dr+ 

+ C||F(«i)-F(u a )|| i - (fl) / -^dr 

Jo r{ u V 

< C + C f R /M , du2/dr dr I C 
A/2 x /2F(u 2 ) 

< C + C (F- 1/2 (l/2) - F~ 1/2 (i?)) < C. 

This proves the lemma for radial solutions. To obtain the estimate in the 
general case, we may always assume that u 2 is the minimal solution of 
JU, so that u 2 < ui and u 2 is radial. By Lemma [2.4l up to replacing / by 
/ given by Remark 12.11 there exists another radial solution V such that 
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V > u\ > U2- Using assumption (|6]), we have 

|/(«i)-/(«a)| dr < [ (/(«i) - /(«a)) dr + JiC||«i-Ua|U-(B) 
(f(V) - f{ua)) dr + K\\ Ul - « 2 || L oo (B) 
< / (/(V)-/(«a)) dr + 2K\\u 1 -u?\\ L ~ {B) 



By Q, ui — U2 is bounded and the result follows from the radial case. □ 

Proof of Lemma l4.2l (|35[) is standard : write the representation formula 
(|34|) both in _B_r and in B\ , change variables in the B\ integral and identify 
the kernels. Next, we prove that given any r £ (0, 1), J gB Gi(x, -)dcr < 1. 
It suffices to show that for any tf> £ C c (0, 1), 

(37) £ <Kr) Qf^ Gi(x,-)doj dr < ||0|| L i( O ,i)- 

By definition of Green's function, the left-hand side of the above inequality 
is the function v solving 

-Aw = <j> in B, 

v — on dB. 

The above equation can also be integrated directly : 
t)'(r) = r 1 ~ N r m^' 1 dt, 



whence \v'\ < \\4>\\l~>-(o,i) an d M < II <^ II (o, i) i- e - lf37)l holds. This proves 
that J gg Gi(x,-)da<l. 

We turn to the second estimate in (|36[1 . Recall that the Green's func- 
tion in the unit ball is expressed for x,y £ B , x ^ y, by 
(38) 

Gi(x,y) = T ((i? 2 + r 2 -2i?rcosv5) 1/2 )-r ((l + R 2 r 2 - 2Rr cos ip) 1/2 ) , 

where R= \x\, r — \y\, ip is the angle formed by the vectors x and y and Y 
is the fundamental solution of the Laplace operator. Differentiating with 
respect to R, we obtain for some Cjv > 0, 

(39) C N ^-(x,y) = 

R — r cos tp Rr 2 — r cos ip 



(R 2 + r 2 ~- 2Rr cos <p) N / 2 (1 + R 2 r 2 - 2Rr cos <p) N / 2 
R — r + r(l — cos op) i?r 2 — r + r(l — cos ip) „ 



((B - r) 2 + 2Br(l - cos<p)) N/2 ((1 - Br) 2 + 2Rr(l - cos<p)) N/2 



cdl9mars2009.tex 



March 20, 2009 



We estimate A and leave the reader perform similar calculations for B. 
Clearly, given e > 0, the expression (|39p remains uniformly bounded in 
the range 1/2 < R, r < 1, e < tp < 2tv — e. Hence, 

\A\ da<C £ +C I \A\da. 
as,. JdB r n[0«p<E] 

For y € 8B r PI [0 < (p < e], let 2 = z(y) denote the intersection of the 
line (Oy) and the hyperplane P passing through x and tangent to the 
hypersphere OBr- Then, there exists constants Ci,C2 > such that for 
all y e dB r n [0 < cp < e], 

Cl(l — COS <j>) < I 2 — x\ 2 <C2(1 — COS (^)). 

Hence, letting B N ~ 1 (x, p) C P denote the N— 1-dimensional ball of radius 
p > centered at i, we obtain 



f \A\ da < C ( 1 + f 

J dB r \ J B- 



\R — r\ + Cr \z — x\ 2 
-l(x,Rsine) (\R - r\' 2 + c\z - x\ 2 ) N ' 2 



lBK-Ho,R e ) (\R~r\ 2 +c\z\ 2 ) N/2 

\R-r\+C\R-rf\£ „_ t 
/ B —(o,^) \R-r\ N (l + c\z\ 2 ) N/2 

I ~ ^—^dz + lR-rl [ \z\ 2 ~ N dz 

l"" 1 (1+C|2| 2 ) 1 is«-l^ ,_Jk_) 

< a 



Working similarly with the B term in (|39[1 . we finally obtain the desired 
estimate (136 D - □ 



5 Asymptotics to all orders 



This section is devoted to the proof of Theorem 1 1.8 1 Our first task consists 
in applying the Fixed Point Theorem to the functional N defined for 
v G B(vo, p), U € I by 



(40) 



[Af(v)](u) = j2 (F(u)-(N-1) 



dt 



where r is given by (|10jl . Let us check first that N (B(vo, p)) C B(vo,p). 
Take v G B(v ,p). Then, 



(41) 



1 > r > 1 - 



— dt = 1 

1 -PJu vo l-p JUQ 



•IF 



dt. 



By ((3]), it follows that for p < 1/4 and £/o sufficiently large, 1 > r > 1/2. 
Hence, 



v , 
- dt 
r 



< C / V2F dt = o{F(u)), 

u 
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where we used Lemma 13,41 So for Uo large and u > Uo, 
N-l 



F(u) 



r. r n 



v , 
- dt 
r 



<P- 



We deduce that 
Af(v) — vo 



(42) 



= 1-4/1 



!'() 



N-l 



Ju o r " 2 F(u) J Uo r 



< P- 



Next, we prove that N is contractive. Given v\,V2 £ B(vo,p), let r\ — 
r(u, Vi),r2 — r{u,V2) (where r is given by ([TUJ ) ■ Then, by estimate (|4ip . 
1/2 < n,T2 < 1 and 



Af(vi)-AT(v2) 



N — 1 



Vl 



dt - 



N - 1 



fa 



< C 



F(w) 
C 





Vl 


V2 




ri 


T2 



< 



F{u) 

C 
F(u) 

C 
F(u) 

Cp 
F(u) 



\vi — Va| dt + 



ru 


1 


1 




\ v 






'u 




^2 





>2F dt + 



I IF dt + 



>2F\ri - r 2 \ dt 



f'i. 



'2-F 



/2_F dt + 



'2F 



'2F 



Vl V2 

ds dt 



dt 



f-U 

J Un 



2F dt. 



F(u)j Uo 

Using Lemma [3.41 we conclude that M is contractive in B(vo, p) if Uo was 
chosen large enough in the first place. We may thus apply the fixed point 
theorem. 

So, it only remains to prove (|13p . We first observe that the sequence 
(«fc) defined by is asymptotic i.e. Vk+i(u) = Vk(u)(l + o(l)), as 

u — » +oo. Since vt + i = J\f(vk), it suffices to prove that Af(vo) — vq = o(«o) 
and iterate. By (@3J, 



A/"(wo) — wo 



vo 



< 



c 

Fin) 



>2F dt 



and the claim follows by Lemma [3.41 So, the sequence (v^) is asymptotic 
and so must be the sequence (ut) defined by (|12|). We are now in a position 
to prove (|13[) . By Theorem 11.11 we may restrict to the case where u is 
radially symmetric. Let v = du/dr. By (|23[) . v solves 



dv 

— H v = f(u) 

dr r 



Use the change of variable u = u(r) to get 
dv N-l 



du 



v = f(u). 
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Integrating, it follows that for some constant C 



C« JV-1 



(43) — = F(u) +C- I -v dt. 

Up to replacing F(u) by F(u) — F(u) + C (which is harmless from the 
point of view of asymptotics) , we may assume C = 0. So it suffices to 
prove that v £ B(vq, p) to conclude that v coincides with the unique fixed 
point of N, whence (|13fl will follow. By (|43[) (with C — 0), v < vq and so 



< v - v < y/2F(u) - ^2 (f(u) - N 1 v dtj 

f" y^Fdt 



< c- 

By Lemma [3.41 it follows that 



Vq — V 

< < v o — v < p 

and v £ B(vo,p) as desired. □ 



6 Universal blow-up rate 



In this section, we prove Theorem 11.91 that is we characterize nonlineari- 
ties for which the blow-up rate is universal. 

Proof of Theorem [TT9] . 

Step 1. We begin by establishing the theorem when Q — B is the 
unit ball. In light of Theorem 11.11 it suffices to prove (|16[1 for one given 
solution u of {l]), which we may therefore assume to be radial. By (|27|) . 
we have after integration that 



(44) f 



1 = dt = 1 - r. 

V2(F(t) - g) 

By definition of uq, we also have 

r + oo -i 

(45) / ; dt = l- r. 

Ju s/2FjFj 

Observe that u > uo, split the integral in (I45p as = + J^°° and 

equate (|44|) and (|45|) . It follows that 



f u 
J un 



dt = — , dt 



^mfi Ju \^2(F(t)-g) ^2F(t) 



f 

■J U 



f 

■J u 



V2F(tj- y/2(F(f)-g) dt 



5 dt 



V / 2F(I)) ^2(F(i) - 5 ) v^W + VWW) ~ 9) 
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Recall that by Lemma[3]4] g — o(F) as t —> +oo and g(u) ~ (N— l)G(w) 
(TV — 1) J Q " df. So, for sufficiently large values of u, we deduce 

(46) L^m*-°f. WWW dL 

Since F is nondecreasing, it follows that 



Hence, 



0<u-uo<CVm^j u J ° 2F(t))3/2 dt, 

and ([in]) follows from (fl3)l . 

Step 2. Next, we prove that (|16[) holds for general domains fi. To this 
end, we combine a standard approximation argument by inner and outer 
spheres (see e.g. [7]) and the comparison technique of [3]. Let u denote a 
solution of (|14[) and take a point xo £ dQ. Let B C Q denote a ball which 
is tangent to dQ. at xo- Shrink B somewhat by letting B £ — (1 — e)-B, 
e > 0. Observe that u G C(B E ) is a subsolution of 



(47) 



At/ = /([/) in 73 £ , 

U — +oo on 9B E , 



By Lemma [231 there exists a solution V E of Q47p . such that V e > u in B E . 
Furthermore, V E can be chosen to be the minimal solution of (|47p such that 
V £ > u in B £ . In particular, V e is radial and e — > 14 is nondecreasing. In 
addition, e — » V e is uniformly bounded on compact subsets of B (working 
as in the proof of Lemma l2.3p . so V £ converges as e — > 0, to a solution V 
of g]) such that V > u in B. By Step 1, 

Jim K(x) - m (1 - rfs(a;)) = 0, 

where ds denotes the distance to dB. Since V > u and since the above 
discussion is valid for any point xo £ dQ, we finally obtain 

(48) Iimsup[u(x) - u (l - d(x))] < 0, 

x— i-dCl 

where d(x) is the distance to dQ. Choose now an exterior ball B C R N \Q 
which is tangent to dQ at xo. For e > small and R > large, the annulus 
A s = i?B \ (1 — e)B contains Q. Let C/ e denote a large solution on A £ , 
which we may assume to be minimal, radial and bounded above on Q 
by u. Again U £ — > U as e — > where [/ is a radial large solution in 
yl = _RB \ B D Q. Repeating the analysis of Step 1. (which was purely 
local) for the case of a radial solution defined on an annulus rather than 
a ball, we easily deduce that 

lim U(x) - u (l - d B {x)) = 0. 

16B 
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Since u> U and since the above discussion is valid for any point xo £ dQ, 
we obtain 



(49) 



liminf [u(x) — Uq(1 — d(x))] > 0. 

x— >dQ 



So, by (|49|) and (|48[) , we have that (|16|) holds in any smoothly bounded 
domain Q. 

Step 3. It only remains to prove that (|16[1 fails when (|17p holds. We 
use Theorem ll.8l to compute the second term in the asymptotic expansion 
of a solution. By (|lip. 



v 1 (u) = J 2 [F(u)-(N-1) / V2F di(l + o(l)) 



- i r° , ; i+oi 



= [ 1 - (JV - 1) ^^1 (1 + o(l)) 



whence 



«i xAFCu) 
1 



1 + (JV-1) 



s (1+0(1)) 



2F( U 



f u v / 2~F lit 



Integrating ()12|) for fc = 1, it follows that for r close enough to 1, 



(50) 



dt 

/2F 



+ (JV-l)(l + o(l)) 



/„' V2F ds 

(2F) 3 / 2 



dt = 1 - r. 



Recall (|45|l . split the integral in (|45|l as J^ 00 = /" an d equate 

(f50l) and |g5j| to get 



dt 

J2F 



= (JV-l)(l + o(l)) 



Since F is nondecreasing, we deduce that 
(51) "iZ^ >(JV-l)(l + o(l)) 

Note also that 

fiV2Fds 



/„' \/2F ds 

(2F)3/2 



+°° /„' V2F ds 
(2F) 3 / 2 



dt. 



dt. 



(52) 



dt< r±dt<^-^ 2 



2F 



Assume by contradiction that lim, (1 - {u\— uo)(r) = 0. Then, (|52|) implies 

that 

r«l V2F ds ( 1„ - 7/n 
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Using this information in (|51[) , we obtain that 



y/2F(u ) Ju (2F) J / 2 

But (|17[) would then lead us to a contradiction with the assumption 
lim,.,^!- (iti — Uo)(r) — 0. So we must have 

liminf (ui — Mo) > 0. 

r->i- 

and so (|16|) fails. □ 



7 The first three singular terms 

In the previous section, we characterized nonlinearities for which only one 
term in the expansion is singular. In the present section, we calculate 
implicitely the next two terms in the expansion. We have not tried to 
characterize those / for which all remaining terms are nonsingular, but 
this can certainly be achieved. We leave the tenacious reader try her/his 
hand at this computational problem. 

We begin by calculating the leading asymptotics of vi, v%. By 
we have 

2 ru 

v J-=F-{N-l) J VS? (l + o(l))d*. 

So, 

1/2 



f 2 =VFll-(N-l)^^(l + o(l)) 



In other words, 



r- N - 1 f u V2Fdt, , , , 



. r V2F dt , 

D— 



Vl = V2F - (N - 1) J (l + o(l)). 

V 2-T 

To calculate «2, we introduce some notation. Given a positive measurable 
function v, set 

f" / ,+ °° dt 

Pv= vdt, Qv= , Rv= — and Tv = (N-l)PQv+P(vRv). 

J v J u v 

v\ is then expressed by 

vi = v -(N-l)(l + o(l))Qv , 

while V2 is given by 

2 /"u 

2 V i + (i)) 



f — (AT — 1) / (u - (JV - l)Qw + o(Q(« ))(l - Rv + o(Rv )) dt = 
= F-(N- l)Pv + (N- l)Tv (l + o(l)). 
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So, 

v 2 = (2F - 2(JV - l)Pv + 2(JV - l)Tv (l + o(l))) 1/2 = 

= «b (l - - + (JV - 1)^(1 + o(l))) = 

jy-l Ppp JV — 1 Tup 3 n2 /^o\' 



+ 



o(TV-F + (iWff ))• 



And so, 

_L = _L N-lPvp N-lTvp 5 

w 2 u I 2 F 2 F 8 



_L_ 

Vo 



+ o{Tv /F+{Pv /F) 2 )^ 
Pvo , ^ ( Tv , 5, 



(Pvo) 5 



: + (N~ 1)- 



3/2 



+ 



(iV-1) 
(2F) 3 / 2 



/2F v ' (2F) 



'2F 



/2F 



+ ■ 



■ (w -i)(r^n 

4 2F v w; 



Integrating once more, we finally obtain 



1 - r 



u 2 (r) \/2~F 



+ (iV-l) 



50 /* V2F dt 
U2{r) (2F)3/2 



du + (1 +o(l))x 



(JV — 1) / - / I (JV 

J«2W V 7 V 



- i 



+ 



'2F 



+ V2F 



5(JV-1) 



(r 



2F 



ds 

72F , 



dt+ 



du 



(2F)3/ 2 ' 



This proves Proposition 11.131 



8 An example: f(u) = u p , p > 1 

Finding the n-th term in the expansion for abitrary n £ N is out of reach 
for general /, simply because of the algorithmic complexity of calculations. 
However, when additional information on / is available, one can guess the 
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general form of the expansion and then try to establish it. This is precisely 
what we do in this section, with the nonlinearity /(it) = u p , p > 1. 

For notational convenience, we shall work with F(u) = \u 2q , where 
2q — 1 = p, which simply amounts to working with a constant multiple of 
the original solution. 

Recall (|lip and (|12p . We want to prove inductively that there exists 
numbers at, bk depending on k,p, N only such that 

n 

(53) v n = u q J2 b k u' kiq - 1] + o(u q - n{q - 1) ), 

k=0 

n 

(54) ? in = d-^^a fe d fc + J B n (d-^ + " +1 ), 

k=0 

where E n (cT~^ T+n+1 ) ~ e n d _ '5 TrT+n+1 for some e n £ R, as d — » + . We 
have uo = V2F = it 9 . Solving for uo in (|12p yields uo = c d~ i-i . So, (|54p 
and (|53|) hold for n = 0. Suppose now the result is true for a given n € N. 
In the computations below, the letter Ck denotes a number depending on 
k,p,N only, which value may change from line to line. By (|53p . we have 

— = u - q ( i + y fefcw-^- 15 + ou-^-v 
\ ti 

= u' q (j2c k u~ k{q - 1) +o(u 
\k=0 



-n(q-l)\ 



So, 

(55) [ + °° ^=t 1 -«fv Cfc t-*(»-^ +0 ( i -("+ 1 )(«- 1 )) 

(56) = ]T erf"*''" 1 ) + (*-(»+ 1 )(«- 1 )). 



It follows that 



1 _ r+°o ds_ 
it v n 



n+1 

_(n + l)( g _l). 



= X)c fc t-* (, - 1) +o(t- 



Whence, 



, _ r-+oo ds_ 
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And so, 



It v n 



u 2q _|_ u q + l ^ Cfc U- fc (9- 1 ) + o(u 1 + 1- n (l- 1 '>) 
fc = 

l+> CfcW + o(l( 

n + 1 



fc=l 

n + 1 



This proves (|53|) . Integrating (|12[) . we obtain 
(57) / — = d. 

Now, w„ +1 =v n + c n+1 u q -^ +1 ^- 1) (1 + o(l)). So, 

— = - + c^-*-^ 1 ^- 13 (1 + 0(1)). 
Vn+1 V n 

It follows that 

/ ,+ °° du f +x du _(„ +2 )(g_i), , 

tf = / =/ — + C„+1M„| 1 ; (l + o(l)). 



In addition, w n ~ «o, so w„ ~ wo, and so U~i j ~ d. Using this in the 
above equation, we get 

d+ Cri+1 d" +2 (i+ (i))= r°°— . 

Recalling that v n is defined by (|57p and satisfies (|54|l by induction hy- 
pothesis, we conclude that 

i n 

«„+! = (rf + C„ + 1 d" + 2 (l + 0(1))"^ Y, ( d + Cn + ld n + 2 (l + o(l)) fc +£ n (d-i i T + " + 1 ). 

fc= 

Expanding again the above expression, we finally obtain 

n + 1 



i)n+i = d i- 1 Y, a kd k + E n+1 (d q-i +n+2 ), 



fe=0 



which proves (|54p . Proposition [T7T4] follows. 

Acknowledgement. The formulation ([2} concerning the positivity of / 
is due to H. Brezis. The authors gratefully thank him for this suggestion, 
and for discussing a preliminary version of this work. 
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